Temperature Resistant Optimal Ratchet Transport 
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Stable periodic structures containing optimal ratchet transport, recently found in the parameter 
space dissipation versus ratchet parameter [PRL 106, 234101 (2011)], are shown to be resistant to 
reasonable temperatures, reinforcing the expectation that they are essential to explain the optimal 
ratchet transport in nature. Critical temperatures for their destruction, valid from the overdamping 
to close to the conservative limits, are obtained numerically and shown to be connected to the 
current efficiency, given here analytically. Results are demonstrated for a discrete ratchet model 
and generalized to the Langevin equation with an additional external oscillating force. 
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Ratchets systems are remarkable profitable due their 
ability to displace particles without an external bias, 
which is known as ratchet transport. One of the chal- 
lenges in the large variety of ratchet phenomena, and 
experiments, is to unveil how large currents of particles 
can be attained and controlled. If such displacement of 
particles can be controlled in distinct ratchet devices, it 
should be of great validity for technological applications. 
In this context we mention some experiments which suc- 
cessfully observed ratchet transport, namely cells to con- 
trol cancer metastasis [1], solids and drops using the Lei- 
denfrost effect [2, 3], micro and nanofluids [4], particles 
in silicon membrane pores [5] , cold atoms [6] , among oth- 
ers. In addition, ratchet models appear as natural can- 
didates to explain directed transport in biology [7], and 
have been realized in a wide range of physical systems [8] . 
In the real-world, ratchet systems always suffer environ- 
mental effects, like temperature, for example, which give 
rise to new current dynamics, and so need to be taken 
into account by any comprehensive approach. Therefore, 
to explain such transport in nature, and for technological 
applications, it is of fundamental priority to understand 
the effect of noise (thermal or not) on the efficiency of 
the ratchet current. 

Recent works [9, 10] have shown that, in the absence of 
temperature, optimal ratchet currents (IZCs) occur along 
organized structures, called isoperiodic stable structures 
(ISSs), in parameter spaces which involve the asymmetry, 
amplitude, phase of the ratchet and dissipation. The ISSs 
should play a fundamental role on the current generation 
in nature, since they bound the optimal 1ZC region, de- 
limit current reversal areas, and explain the remarkable 
diversity of possible IZCs. Results have partially been 
extended to quantum systems [11]. The ISSs are Lya- 
punov stable islands and are supposed to be generic in 
dynamical systems [12-15]. Thus, any attempt to un- 
derstand general effects of the temperature on the 7ZC, 
should be performed by analyzing the behavior of the 



ISSs, and their surroundings, when the temperature is 
taken into account. This is the main goal of the present 
Letter: to show the transformations of the ISSs under 
temperature, and that optimal IZCs inside the ISSs are 
resistant to reasonable temperatures. This is the opposite 
of what is done in non-inertial thermal rocket ratchets, 
where IZCs are thermally activated. Our results also con- 
firm that generic ISSs, like Shrimps discovered in deter- 
ministic dynamical systems [12, 13], for example, persist 
to some noise intensity in stochastic systems. Up to our 
knowledge, there is only one work [16] which observed the 
tendency that currents (not IZCs) survive to reasonable 
noise effects inside ISSs. 

A model which presents all essential features regarding 
unbiased 7£Cs, is the dissipative map [6] 



Vn+i = 7Pn + K[sin(x n ) + asin(2x n + </>)] + f (n). 



(i) 



where p n is the momentum variable and x n de position, 
n = 1,2, ...,7V represents the discrete time, K is the 
nonlinearity parameter, and £(n) is the stochastic vari- 
able obeying (f (n)) = and (£(n) 2 ) = 2(1 - <y)k B T, 
with ks = 1 being the Boltzmann constant and T the 
temperature. For T = this model was considered for 
K = 6.5, a = 0.5,0 < 7 < 1 in [17], and in the param- 
eter spaces combining pairwise all parameters 7, K, a, <fi 
[9, 10]. The dissipation parameter 7 reaches the over- 
damping limit for 7 = and the conservative limit for 
7 = 1. The ratchet effect appears due to the spatial 
asymmetry, which occurs with a 7^ and (j) ^ mir 
(m = 1,2,...), in addition to the time reversal asym- 
metry for 7 7^ 1. The above model is quite general to 
understand IZCs. The parameter K represents the inten- 
sity of an external asymmetric force in space. But since 
for 7 = 1 the above model can be derived from a kicked 
Hamiltonian, K is as well the external kick which drives 
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the system out of equilibrium, which is a necessary con- 
dition to obtain ratchet transport. 
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Figure 1. (Color online) The 1ZC (see color bar) plotted in 
the parameter space (if, 7) with a grid of 10 3 x 10 3 points, 
10 4 ICs with (po) = (xo) = inside 



a 0.5, tt/2, M 
the unit cell (— 27r,27r), 



N = 10 b ,n c = 9 x 10 4 iterations for 



(a) T = 10" 5 , (b) T = 10" 3 and (c) T = 10" 



For a small temperature T = 10 5 , Fig. 1(a) shows the 



j_ 



E 



M 



iV- 



(colors) as a function 



of 7 and if. M is the number of initial conditions (ICs), 
N is the total number of iterations used to calculate the 
7£C, and n c is obtained empirically and represents a cut- 
off value to avoid transient effects. The black color is 
related to small currents; dark green, green to white col- 
ors are related to increasing positive currents, while red 
to yellow colors related to increasing negative currents. 
Three main regions with distinct behaviors can be iden- 
tified (for more details see Fig. 1 in [9] for T = 0): (i) 
a large "cloudy" background A, mixed with black, dark 
green and red colors, showing a mixture of zero, small 
positive and negative currents. This is a chaotic region 
[10] where the 1ZC is not that large, neither efficient, and 
is independent of the magnitude of the positive Lyapunov 
exponent; (ii) structures Bl, Cl and Dl with sharp bor- 
ders (the ISSs) and distinct colors, which are embedded 
in the cloudy background region (L is an integer or ra- 
tional number). These ISSs are responsible for the larger 
IZCs and they align themselves along preferred direction 
in the parameter space (see [9, 10]). Inside the ISSs the 
Lyapunov exponents (LEs) are zero and negative, and 
the motion is periodic, while a sequence of doubling bi- 
furcations cascade occurs in a specific direction. In the 



Cl and Dl ISSs, these bifurcations start from the in- 
ner part [10] of the ISSs, and reach the chaotic motion 
at specific borders. In the Bl ISSs, saddle-node bifur- 
cations of period- 1 occur on the left and right borders, 
and were determined analitically in [9, 10]. Due to a 
crisis the analytical left borders do not match exactly 
with the numerical ones. A sequence of periodic dou- 
bling bifurcations occur when going down (smaller 7) to 
the right, where K values increase [see white arrows in 
Fig. 1(b)]. It was also shown [9, 10] that the 1ZC is in- 
dependent on the period of the orbits. The values 2tt L 
determine the magnitude and direction of the 1ZCs. For 
example, the 1ZC inside the ISSs Bl, known as cusp- 
idal singularities (cusps), increases positively along the 
sequence B\,B2, . . The 1ZC inside the Cl ISSs, in- 
creases negatively along the sequence C_i, C_2, — The 
smaller Dl ISSs are the shrimp-shaped, well known to 
appear in many dynamical systems. In general, it was 
shown for T = [10] that the 1ZC is always more efficient 
inside the ISSs; (Hi) the last region has strong positive 
and negative currents (region E), with not well defined 
borders, and occurs close to the conservative limit 7 = 1. 
This is related to a chaotic region (as shown in [6] for 
T = 0), where the large 7ZCs are generated due to the 
accelerator modes which exist in the conservative limit. 
In Fig. 1(a) the temperature T = 10 -5 is small and all 
ISSs are mainly preserved. However, the ISSs have small 
elongated antennae, as seen in Fig. 1 from [9], which al- 
ready disappeared for this small T. Consequently, the 
lines connecting the distinct and D_i/ 2 ISSs, were 

destroyed by noise. LEs inside these antennae become 
positive, meaning that the ISSs antennae are destroyed. 

Figure 1(b) shows the same parameter space, but 
now for a temperature T = 10 -3 . Comparing with 
Fig. 1(a), we observe that the lower thin antennae from 
the B2,Bs, . . . ISSs were destroyed. The Cl ISSs begin 
to be destroyed from their borders, which are not sharply 
defined anymore. The Dl shrimp-shaped ISSs have al- 
most disappeared, and the optimal IZCs inside them are 
substituted by the smaller IZCs from the chaotic region. 
The magnitude of the IZCs inside the remaining ISSs is 
essentially not affected by T, and still increases as the 
conservative limit is approached. All missing ISSs, or an- 
tennae, gave place to chaotic regions with smaller IZCs. 
In the chaotic regions A and E, the 1ZC is almost unaf- 
fected by T. 

As the temperature increases to T = 10 -2 , a signif- 
icant amount of ISSs start to be destroyed. This is 
shown in Fig. 1(c). Almost all antennae from the Bl 
ISSs are destroyed, remaining, in general, only the inner 
part of these structures, where the lowest period of the 
ISS is found. The less affected ISS is B±, which originally 
had the smallest 1ZC from the sequence B\ , B2 , — The 
boundaries of these ISSs are not sharp anymore, but ap- 
pear blurred. Remarkable is that the magnitude of the 
7£Cs, inside the inner part of the remaining Bl ISSs, is 
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almost the same. Along the period doubling bifurcations 
[see white arrows in Fig. 1(b)], T destroys the optimal 
7£C, remaining just a region (now chaotic) with IZCs 
values a tiny amount larger when compared to the sur- 
rounding background chaotic region. Besides that, the 
Cl ISSs are completely destroyed, giving place to the 
chaotic motion with almost zero current. Comparing to 
Fig. 1(a), in the chaotic regions A and E the TZCs remain 
unchanged. Close to if ~ 2.5 and 7 ~ 0.9, we observe in 
Fig. 1(c) a region where larger 1ZCs occur. These larger 
IZCs were thermally activated, since they do not show up 
in Figs. l(a)-(b). It gives an example, in distinction to 
all other regions of the explored parameter space, where 
the 1ZC is activated by T, instead of being destroyed. 

In this context the key question is: what is the critical 
T to destroy the 1ZC in different stable regions (ISSs) of 
the parameter space? To answer this question the fol- 
lowing simulation was performed. For each parameter 
combination (if, 7), T was increased until the 1ZC was 
close to zero, within the precision A ~ ±3.0, which is a 
value bellow the minimum 1ZC allowed by the periodic 
constraints in ISSs with main periods 1 and 2 [9], that is, 
the main part of our parameter space. It is possible to 
observe that the close to zero 7ZC in the chaotic region, 
is independent of the parameters (if, 7) within the pre- 
cision A. This defines the critical temperature T c , where 
the optimal 7£C, due to periodic motion (ISSs), is trans- 
formed into the small 7ZC due to the chaotic motion. 
Figure 2(a) shows T c (color bar) in the parameter space 



Figure 2(b) shows in details T c from the box of Fig. 2(a), 
with three connected shrimp-like structures. Thus, the 
ISSs can still be recognized in such T c analysis. Even 
though T c is larger inside the inner parts of the ISSs, it 
tends to decrease at their borders. 

The physical origin of the thermally induced destruc- 
tion of the optimal 1ZC can be rather complex. It depends 
very much on which region in the parameter space is con- 
sidered. Inside the ISSs there is always at least one pe- 
riodic stable attractor, which coexists with the unstable 
periodic orbits and chaotic repellers. With the introduc- 
tion of T, the attracting motion (periodic or not) becomes 
transient and the dynamics is ergodic, as earlier studies 
about noise effects and transient chaos in dynamical sys- 
tems indicate [18]. Temperature also induces transitions 
between all possible dynamics: periodic/non-periodic at- 
tractors, chaotic repellers etc. Thus, for T ^ 0, the final 
IZCs will be an average over the velocities of each of such 
dynamics, weighted by their lifetimes. Since the 1ZC is 
more efficient inside the ISSs [10], we should understand 
the effect of T on the 7ZCs efficiency, which is defined by 
V = ((P)) 2 /\((P 2 )) ~ A)| [19]. D = k B T is a measure of 
the thermal energy and ((•••)) denotes an average over 
ICs, thermal realizations and time. For the ergodic case 
it reduces to just one average, over the asymptotic dis- 
tribution. Considering that the thermal "kicks" are typ- 
ically weak enough, so they can be neglected after one 
iteration, but rarely strong enough to drive the orbit out 
of the attractor, then it is possible to show that 




Figure 2. (Color online) The critical temperature T c (see color 
bar) necessary to destroy the optimal 1ZC plotted in the pa- 
rameter space (a) from Fig. 1 and (b) inside the box from 
Fig. 2(a) and for a = 0.5, <j) — tt /2 and same ICs from Fig. 1. 



In (c) we show that T c 
7c 0.5. 



|7 — 7c I , for if — 2tt and around 



(if, 7). T c increases from black (T c ~ 0.0) to red, yellow 
and white (T c ~ 0.07). It is clearly observed that T c is 
larger inside the ISSs, especially inside their inner parts. 



£^<P?r=o> + (! " 7)A)/(1 + 7) 



(2) 



where the sum is over all attr actors, (p^ T=0 ) 1S the 
mean of the momentum square on the attractor i (for 
T = 0), OLi = \i{Ti is the statistical weight of each at- 
tractor, where T{ is the mean lifetime on attractor z, and 
Hi = linin^oo (rrii)/n, where mi is the number of times 
a given trajectory visited the attractor i during the time 
n. According to our simulations, transitions between 
deterministic dynamics [the sums over attractors from 
Eq. (2)], induced by T, affect more the efficiency than 
pure thermal effects (proportional to Do). Such tran- 
sitions may even increase the efficiency by allowing the 
orbit to access attractors with higher currents related to 
the unstable dynamics from T = 0, as could be seen in 
Fig. 1(c). However, inside the ISSs the current is usually 
optimal, thus T tends to decrease its efficiency. Also, as 
T c is related to very frequent transitions between differ- 
ent dynamics, the efficiency for larger T must decrease 
substantially. 

As far as we have knowledge, there is no first principles 
theory that quantifies the critical temperatures which 
destroy the 7£Cs. Quasipotentials [18] and escape rates 
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[16] from periodic/non-periodic attractors and chaotic re- 
pellers, have been used to quantify T c in a context not re- 
lated to 7ZCs. It is known [18] that close to a saddle-node 
bifurcation, the dynamics of a higher-dimensional system 
can be reduced to an one-dimensional normal form, and 
then the critical Gaussian noise follows a c ~ |7 — 7 C | 3 ^ 4 , 
with <j c being the critical noise amplitude. To check this, 
we plot T c in Fig. 2(c) as a function of I7— j c \, for K = 27r 
and 7 C = 0.5. This is a point at the border line from the 
B\ ISS (see white line in Fig. 2(a)), which was calculated 
analytically in [9, 10], and has a saddle- node bifurcation. 
Since in our case a oc y/T we expect that T c ~ I7 — 7c| 3 ^ 2 - 
This is exactly what is observed in Fig. 2(c) from our nu- 
merical simulation (cross points). For comparison, we 
plotted the exact continuous line |7 — 7 C | 3 ^ 2 - 

To generalize our results to a continuous, and more 
realistic system, we analyze the T c to destroy the 
optimal 1ZC in the Langevin equation: x + jx — 
5.0 [sin (x) +0.7 cos (2;r)]-if t sin (£)+£(£) = 0, where K t 
is the amplitude of the external time oscillating force, 7 is 
the viscosity, while the ratchet potential and the stochas- 
tic term ^(t) are like from Eq. (1). ISSs have been found 
[9] for this problem for T = in the parameter space 
[K t ,x — e^ -7 ^]. Figure 3 shows the critical temperature 
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temperatures. The ISSs persist in such a plot. Distinct 
regions in the parameter space present different dynam- 
ics and T affects the optimal IZCs distinguishable. The 
general physical explanation is that T induces a tran- 
sient chaos and also transitions between all possible dy- 
namics: periodic and non-periodic attractors, chaotic re- 
pellers etc. An analytical expression for the 1ZC efficiency 
was obtained, which combines all these dynamics. It es- 
sentially shows that there are two thermal effects on the 
efficiency, namely the slight variation in (p 2 ) inside the 
attr actor, and the transitions between different dynam- 
ics. Considering the scale of T c , the first effect is practi- 
cally neglectable and thus the 1ZC is optimal inside the 
ISSs, where the transitions are very rare. For T^0, the 
periodic attractor loses its dominance on the 7ZC, allow- 
ing the rise of those new kinds of dynamics that control 
the 1ZC above a critical value of T, T c . This tempera- 
ture is obtained numerically in Fig. 2 for the asymmet- 
ric standard map, in Fig. 3 for the Langevin equation 
with an external driving force, and is a direct measure 
of all complicated combined effects to destroy the 1ZC. 
It also gives quantitative information about the "global 
stability" of attractors under noise in phase space. This 
allows us to relate the robustness of the optimal RC to 
the global stability of solutions, rather than to the local 
stability obtained from the Lyapunov analysis. As well 
motivated by [20], information regarding global stability 
is usually more relevant in real- world systems where per- 
turbations are finite. Hence, our results underscore the 
fundamental role of isoperiodic stable structures in the 
real- world optimal ratchet transport, where environmen- 
tal effects cannot be neglected. 



Figure 3. (Color online) The critical temperature T c in the 



parameter space (K t ,x 



B (-7h 



of the Langevin equation. 



The shrimp-like ISS can be identified. 

T c , required to destroy the 1ZC in a portion of the pa- 
rameter space (Kt,x). The black color (T c = 0) means 
that no optimal 1ZC is observed, red, yellow to white col- 
ors are related to the increasing values of T c required to 
destroy the 1ZC. The shrimp-like ISS can be identified, 
showing that its inner part is more resistant (yellow) to 
external noise, while on its border (red), smaller noises 
are enough to destroy the optimal TZCs region. 

Concluding, optimal ISSs IZCs in the parameter space 
are shown to be resistant to reasonable temperatures 
T, i.e., it exists for thermal energies around 1000 times 
smaller then the transport energy. Increasing the values 
of T, the IZCs inside the ISSs begin to be destroyed from 
their borders, and the ISSs start to become more and 
more blurred. We determined numerically the critical 
temperature T c , required to destroy the optimal 1ZCs in 
the whole parameter space. Results show that, in gen- 
eral, as nearer of the ISSs' center, higher are the critical 
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